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ABSTRACT 

Observations show that the solar p-mode frequencies change with the solar 
cycle. The horizontal-phase-velocity dependence of the relative frequency change, 
scaled by mode mass, provides depth information on the perturbation in the solar 
interior. We find that the smoothed scaled relative frequency change varies along 
| the solar cycle for horizontal phase velocities higher than a critical value, which 

j~ J | corresponds to a depth near the base of the convection zone. This phenomenon 

suggests that the physical conditions in a region near the base of the convection 
zone change with the solar cycle. 
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1. Introduction 
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Observations show that the solar p-mode frequencies vary with the solar cycle, and 
the change in frequency correlates with magnetic activities. The previous study suggests 
that the measured frequency change is caused by the change in physical conditions near 
the surface (Libbrecht & Woodard 1990; Shibahashi 1991), and no evidence of solar-cycle 
variations of physical conditions deep in the solar interior has been found from the measured 
frequencies (Howe et al. 1999; Basu & Schou 2000; Antia et al. 2001; Basu & Antia 2002; 
Eff-Darwich et al. 2002; Basu et al. 2003). It is expected that the perturbation in the solar 
interior have a rather small contribution to the frequency change, if it exists, because of 
small 1/(3. It is also reasonable to expect, based on the current dynamo theories, that the 
second largest contribution to the frequency change is from the perturbation near the base 
of the convection zone (BCZ). A perturbation near the BCZ would modify the frequencies 
of different modes in different ways: the frequencies of the modes which can penetrate 
into the BCZ would change, while the modes which can not penetrate into the BCZ would 
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remain intact. The penetration depth (lower turning point) of a mode is determined by 
the horizontal angular phase velocity w (= cu/L = + 1)] 1//2 ), where to is the angular 

frequency and / the mode degree. Thus the contribution of the BCZ perturbation to the 
frequency change depends on w. Here we present the evidence that the w dependence of 
smoothed scaled relative frequency change, (5u> n i/uj n i)E n i, varies with the solar cycle at large 
w, while remains constant at small w, where 5u) n i is the frequency change and E n i the mode 
mass. 

To search for the contribution of the near-BCZ perturbation to the frequency change, 
we have to remove the contributions of the near-surface perturbation, which is much greater 
than that of the BCZ perturbation. In this study we adopt a method similar to one de- 
veloped by Di Mauro et al. (2002) to remove the near-surface contributions to solar-cycle 
variations of frequencies. This method was originally developed to remove the contributions 
of uncertainties near the surface in solar models to the p-mode frequencies. Here we use it 
to remove the contributions of the near-surface perturbation caused by solar activities. 



2. Method 

In the asymptotic theory, the frequency change 5u n i, due to sound-speed perturbations 
can be expressed as (Di Mauro et al. 2002) 

Six) 

— —Snl = Hi(w) + H 2 (u n i) + H 3 (u n i, w) + H A (u nh w) (1) 

where S n i is approximately proportional to the mode mass E nt . In this study we replace S n i 
by E n i = Eni/Ezifi, where £21,0 is the mode mass for mode n — 21 and 1 = 0. In equation 
(1), Hi is from the perturbations in the whole Sun, while H 2 , H 3 , and if 4 , are caused only by 
the near-surface perturbation. Among the near-surface terms, H 2 (ui n i), is the lowest order 
term. It depends on frequency only, and could be expressed as (Di Mauro et al. 2002) 

M 

H 2 (u nl )=J2 a ? )p ^ ( 2 ) 

i=i 

where Pi are the Legendre polynomials, and u n i is the frequency scaled to be in the interval 
[— 1, 1] for the modes considered. The third and fourth terms, H 3 and H±, are one order 
higher than H 2 (Brodsky & Vorontsov 1993; Gough & Vorontsov 1995). They are functions 
of uj n i and w, and can be expressed as (Di Mauro et al. 2002) 

1 M 

H 3 (u nl ,w) = — VV {3) Pj(u; n ;) 
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The first term in equation (1), Hi, is from the perturbations in the whole Sun. It is a 
function of w only, and can be expressed as an integral (Christensen-Dalsgaard et al. 1988) 
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where r t is determined by r t = c(r t )/w. The largest contribution to Hi is from the near- 
surface perturbation. The second largest contribution is expected from the perturbation 
near the BCZ. Here we define a critical horizontal phase velocity Wb such that the modes 
with w < Wb are not affected by the BCZ perturbation, while the modes with w > Wb are 
affected by the BCZ perturbation in the ray theory. It is reasonable to assume that the 
only non-negligible contribution to H 1 for the modes with w < Wb is from the near-surface 
perturbation. If the near-surface perturbation is located within a thin layer at r , from 
equation (4) Hi can be approximated by 



Hi(w) pa s 



w 2\ -i/2 
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w z 



for wo < w < Wb 
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where Wo = c(r )/r , and s = (Sc/c 2 ) dr. The perturbation near the BCZ has a contri- 
bution to (5uj n i/uj n i)E n i only for the range of w > Wb through Hi. In the following analysis, 
we limit our study in a range of 2.5 < ujjl-n < 3.5 mHz and 190 < w < 1570, where w is 
computed with frequency in units of /zHz. For this range, (wo/w) 2 < 10~ 2 , if ro > O.993i? . 
Thus Hi(w) from the near-surface perturbation can be further approximated by 



Hi(w) pa so 1 + 
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The contributions from the near-surface perturbation can be expressed as 
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If we define [(5u n i/u; n i)E n i]i, cz to be the contribution of the perturbations near and below 
the BCZ, (5uj n i/uj n i)E n i can be expressed as 
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where [(5u n i/u n i)E n i] hcz = for w < w b . 

To obtain [(Su n i/u n i)E n i] bcz , first, we fit equations (7) to the measured {5uj n i/uj n i)E n i for 
the range of w < w b to determine the coefficients, af\ af \ af \ and af\ Second, we use 
these coefficients to compute [(5u> n i/uj n i)E n i} SVLr f for all ranges of w. Third, the contribution 
from the BCZ, [(5uj n i / 'u n i)E n i}^ cz , is computed with equation (8). 

3. Modeling 

In this section, we test the above method with a model computation. First, we perturb 
a standard model (Christensen-Dalsgaard et al. 1996) by changing its adiabatic exponent 
Ti(= d\np/d\n p\ s ) in two regions: one is located near the surface and another near the 
BCZ. Each perturbation is a Gaussian in depth. The near-surface perturbed region has a 
magnitude of 5T 1 /T 1 = 5.95 x 1CT 3 located at O.9995i? with a width of lCr 3 i? . The 
magnitude of STi/Ti is chosen such that the frequency change is close to the measured value 
at maximum. The near-BCZ perturbed region has a mag nitude of 5iy r i = 4 x 1CT 5 located 
at O.713i? with a width of O.O5_R . Second, the mode frequencies of the perturbed model are 
computed (Christensen-Dalsgaard 2003). Third, the relative frequency difference between 
the perturbed and unperturbed models scaled by mode mass, (5uj n i/u) n i)E nh is computed. 
Its value versus log(w) is shown in the upper left panel of Figure 1. The scattering of the 
data is due to the frequency dependence. It is noted that an increase in {5uj n i/uj n {)E n i is 
visible around log(w) = 2.7 which corresponds to the location of the assumed perturbation 
near the BCZ. To see the trend of the scattered data, we smooth data with a 41-point box 
running mean. The advantage of smoothing is discussed in §5. The oscillatory features in the 
smoothed data are due to the ridge structure of p-mode frequencies. Fourth, equation (7) is 
fitted to (5u n i/uj n i)E n i for the range of w < w b to determine the coefficients af\ af \ af \ and 
af\ Here we adopt w b = 430 (log(w) = 2.63), and M — 8. Fifth, these coefficients are used 
to compute [(Suj nl /u n i)E n i} SUT{ with equation (7) for all ranges of w. Sixth, [(5u n i/u nl )E n i] hcz 
is computed with equation (8). The result is shown in the upper right panel of Figure 1. 

The scattering of data is greatly reduced, especially for w < w b . It is because the fre- 
quency dependence, described by [(5u n i/u n i)E n i] SUT f, has been removed. The global trend of 
increase in {8uj n i/uj n i)E n i with w is also removed. Its w dependence is flat and its value is close 
to zero. For comparison, we plot (5u> n i/u n i)E n i computed with only the BCZ perturbation 
in the lower part of the upper right panel in Figure 1. It can be seen that the near-surface 
contributions are reasonably well removed and a small change in frequency caused by the 
near-BCZ perturbation is reasonably well recovered. 
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4. Data Analysis and Results 

The solar p-mode frequencies used in this study are derived from the data taken with the 
Michelson Doppler Imager (MDI) on board the SOHO (Scherrer et al. 1995). The frequencies 
are measured with a time series of 72 days. We use the modes in the range: 2.5 < u/2n < 3.5 
mHz and 190 < w < 1570, because of no missing frequency in this range. To avoid that the 
result may be dominated by a few frequencies with large errors, we discard those frequencies 
by the following procedure. First, for each mode, temporal variations of frequency is fitted 
to a linear function. Second, the standard deviation is computed. Third, the frequencies 
deviated from the fit greater than three standard deviation are discarded. Less than 1% of 
frequencies considered are discarded. 

The frequency averaged over a solar minimum period, 1996.05 - 1997.07, is used as the 
reference frequency. The frequency change along the solar cycle is computed relative to the 
reference frequency. The relative frequency change, scaled by mode mass, averaged over a 
maximum period (2000) is shown in the lower left panels of Figure 1. The data are rather 
scattered. Part of scattering is due to the frequency dependence of (5u n i/u) n i)E n i as shown in 
the model computation in §2, part is due to the error of measured frequency. The smoothed 
data with a 41-point box running mean as §3 are shown with a solid line. The interesting 
phenomenon shown in the smoothed data is that the values of {8uj n i/uj n i)E n i at small w and 
large w are different. It is approximately constant for logiw) < 2.7 and drops to a smaller 
value at log(w) 2.7. 

We use the method described in §2 and §3 to remove the contributions of the near- 
surface perturbation to obtain [(Su n i/uj n i)E n i] hcz . The results are shown in the lower right 
panels of Figure 1. Although the scattering of data is reduced, it is still significantly greater 
than that in modeling. It is probably caused by the the error of measured frequency. The 
smoothed [(5u n i/u n i)E n i]i, cz is also shown by a solid line. In the range of w < Wb, similar 
to the modeling in §3, the global trend in (5uj n i/uj n i)E n i is removed. The w dependence is 
rather flat and the value is about zero. The difference between small w and large w becomes 
more apparent in [{5u nl /u nl )E nl \ hc% . 

To see how this phenomenon varies with the solar cycle, we plot the smoothed (5u n i/uj n i)E n i 
for different periods along the solar cycle in the upper left panels of Figure 2. In the period of 
low activity (1997-1998), the curve is rather flat. As the activity increases, (5u n i/uj n i)E n i at 
log{w) > 2.7 becomes smaller relative to that at logiw) < 2.7. The difference increases with 
the activity. After subtracting the near-surface contributions, [(5cu n i/uj n i)E n i\i, cz is shown in 
the upper right panel of Figure 2. At logiw) < 2.7, [(5ui n i/ui n i)E n i\b cz is approximately flat 
and close to zero. At logiw) > 2.7, [(Scu n i/ui n i)E n i\ bcz is negative. 
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5. Discussion 

The w dependence of (5u n i/uj n i)E n i provides depth information on the perturbation in 
the solar interior. However, it is difficult to see the weak w dependence of (5u n i/u n i)E n i from 
the scattered {5uj n i/uj n i)E n i. We smooth {5uj n i/uj n i)E n i to search for its weak w dependence. 
Smoothing has another advantage. Smoothing is equivalent to averaging over frequencies at 
a fixed w. If the higher order terms, H 3 and H4, are neglected in equation (1), smoothed 
(5u n i/uj n i)E n i equals to H\{w) + constant, where the constant is the average of H 2 over 
frequencies which is independent of w. Thus smoothed {8uj n i/uj n i)E n i gives the w dependence 
of Hi. This can be demonstrated with model computations. 

We have tried several ways of smoothing: the box running mean, fit to a polynomial, and 
Savitzky-Golay smoothing filters (Press et al. 1999). Although different smoothing methods 
yield different results, the main feature discussed here remains similar. Here we use a 41- 
point box running mean to smooth (5uj n i/uj n i)E n i. The w dependence of the smoothed 
(5uj n i/uj n i)E n i varies with the solar cycle as shown in Figure 2. The smoothed (5uj n i/uj n i)E n i 
remains rather flat at low w in all periods, but it begins to change (decrease) at log{w) ~ 
2.7. The difference between small w and large w increases with the magnetic activities 
on the surface. The error bars shown in Figures 1 and 2 are computed with the errors of 
measured frequencies, based on the error propagation. It is noted that the fluctuation of 
smoothed {8uj n i/uj n i)E n i in the range of logiw) < 2.7 is approximately the same for all periods 
and smaller than the temporal variations at log(w) > 2.7. Together with the small error 
bars of smoothed (5cu n i/u! n i)E n i, it suggests that the random errors of measured frequency 
can not account for the temporal variations at log(w) > 2.7. We can not rule out the 
possibility of systematic errors, though it is not clear what kind of systematic error would 
cause this phenomenon. After removing the near-surface contributions, the relative decrease 
at logiw) > 2.7 becomes more apparent. 

Negative [(5u n i/uj n i)E n i]^ cz at logiw) > 2.7 correspond to the existence of a region near 
the BCZ with a smaller wave speed relative to solar minimum. From model computations, a 
decrease in sound-speed near the BCZ would create a dip in [(5cu n i/cu n i)E n i] hcz at log(w) ~ 
2.7. Because of the small S/N, it is difficult to discuss the existence of this dip in Figure 2. 
However, we can increase the S/N by computing the difference between [{5oj n i/u) n i)E n i}\> cz 
averaged over small w and large w. Figure 3 shows the difference between [{Su n i / u n i) E n i\ hcz 
averaged over log{w) = 2.4 — 2.65 and log{w) = 2.7 — 3.0 versus time. This difference 
approximately increases with solar activity. The cause of the relative decrease in frequency 
at large w is not clear. If magnetic fields at the BCZ causes an increase in wave speed, the 
phenomenon observed here might suggest that variations of magnetic flux near the BCZ and 
surface magnetic flux are not in phase. 
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In this study, we use the frequency at solar minimum as the reference and investigate 
solar-cycle variations of frequency change relative to this reference frequency. This approach 
suffers from the risk that the error in frequency at minimum affects the results in all other 
periods. To avoid it, we also use the frequency averaged over all periods as the reference. 
The smoothed (5cu n i/ui n i)E n i, are shown in the lower left panel of Figure 2. Although these 
curves are smoother than those using the frequency averaged over minimum as the reference, 
the relative temporal variations of two groups of curves are consistent. 

Since S n i is proportional to the mode mass E n i, in this study we replace S n i by E n i/E2i i o- 
Although sometimes Q n i, which equal E n i normalized by mode mass of a radial mode of the 
same frequency, are used in other studies, the frequency dependence of the mode mass of 
radial modes would make the first term in equation (1) no longer a function of w only. 
However, in practice we find that using E n i and Q n i yield very close results. The difference 
between [{5uj n i/uj n i)Qni\hc7, averaged over small w and large w using Q n \ is also shown in 
Figure 3. 

In removing the near-surface contribution, the choice of the range of w in the fit would 
affect the detailed features of [(5u n i/u n [)E n i} hcz . However, the main feature of smoothed 
[{8ui n i/uj n i)E n i\\, cz is robust with respect to the choice of the fit range. In deriving equation 
(6) from equation (4), we have assumed that the perturbation in the CZ is negligible except 
near the surface. If it is not the case, we need to include higher-order terms of w~ 2 in 
equation (6). 

In summary, Figures 1-2 indicate that solar-cycle variations of the w dependence of 
(5u n i/uj n i)E n i at log{w) < 2.7 and log{w) > 2.7 are different. This phenomenon suggests 
that the physical conditions in a region near the BCZ change with the solar cycle. If we 
attribute the difference in [(5u n i/uj n i)E n i]i, cz between log(w) < 2.7 and log(w) > 2.7 to 
a change in Ti near the BCZ, it corresponds to — 5r 1 /r 1 (= —5c 2 /c 2 ) ~ 2 — 6 x 10~ 5 at 
r 0.65 — O.67i? , if the FWHM of the Gaussian perturbation is 0.05-R©. 
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Fig. 1.— Top: (5u nl /u nl )E nl (left panel) and [{8u n i/u nl )E nl \ hcz (right panel), computed with 
models, versus log(w). The solid lines are the smoothed data using a 41-point box running 
mean. For comparison, (5u n i/uj n i)E n i computed with only the BCZ perturbation is shown in 
the lower part of the right panel. Bottom: {5uj n i/uj n i)E n i (left panel) and [(5ui n i/ui n i)E n i\i, cz 
(right panel), difference between year 2000 and solar minimum, versus log(w). The error 
bars at log{w) = 2.45, 2.64, and 2.81 are shown. 
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Fig. 2. — Top: Smoothed (5uj n i/uj n i)E nl (left panel) and [(5u n i/u n i)E n i\ hcz (right panel) 
versus log(w) for various periods along the solar cycle, using the frequency averaged over 
solar minimum as the reference. Bottom: Smoothed {8uj n i/uj n i)E n i versus log(w) for various 
periods along the solar cycle, using the frequency averaged over all periods as the reference. 
The period of 1996-1997 is 1996.05 - 1997.05. The period of 1997-1998 is 1997.07 - 1998.10. 
The curves of smoothed (5u n i/uj n i)E n i are shifted to avoid interference: 1996-1997 is shifted 
by 4.5 x 10" 5 , 1997-1998 by 2.5 x 10~ 5 , 1999 by -1 x 10~ 5 , 2000 by -2 x 10" 5 , and 2001 
by —1 x 10~ 5 in both top and bottom plots. The error bars at log(w) = 2.45, 2.64, and 2.81 
are shown. 
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Fig. 3. — Difference between averages over small w and large w versus time for (5uj n i/uj n i)E n i 
(filled circle) and [{Suj n i/uj n i)E n i]^ cz (open circle). The triangles denote the result for 
l(Su n i/u n i)Q n i]b cz . It is noted that there are only seven-month data in 1998. 



